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Abstract 

We develop an Over Barrier Model for computing charge exchange between ions and one-active- 
electron atoms at low impact energies. The main feature of the model is the treatment of the 
barrier crossing process by the electron within a simplified quantum mechanical formulation which 
^ i takes into account: (a) the probability of electron reflection even for over-barrier scattering, and (b) 

_ ^ ' the discreteness of the receiving atom's quantum levels which strongly suppress captures far from 

'c/2 ' 
J>^' the resonance condition. It is shown that inclusion of these effects yields a fairly good prediction 

, ^ , of experimental data. We discuss also the probability of electron re-capture by the target. 
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I. INTRODUCTION 

Charge exchange processes between slow atomic particles are of great importance in 
plasma physics and astrophysics. By "slow" we mean that the interparticle velocity is 
smaller than the classical velocity of the exchanged electron. 

While only quantum mechanical (QM) methods can give really accurate computations of all 
of the basic quantities for these processes, i.e. total and partial or differential cross sections, 
less precise but simpler methods can still be highly valuable when only moderate accuracy 
is sought. In the medium-to-high impact velocity range most preferences go to the Classical 
Trajectory Monte carlo (CTMC) method, which is also more and more often successfully 
applied also to the low velocity range (see e.g. JI], |] for a discussion and some recent im- 
provements on this subject). However, the CTMC method has two disadvantages: (i) it 
is entirely numerical in character, thus somewhat masking the underlying physics; (ii) it 
relies on large numbers of simulations, thus being rather time-consuming. For these reasons, 
analytical or semi-analytical methods can still be useful. Over barrier models (OBM) are an 
example of these models. They are known since a long time and are still being improved 
to include as much physics as possible 0, |^, |^, |^, |^ . 

In this work we present a new version of the OBM. It is based upon the papers @, 0! if) 
but is more than a simple refinement: we adopt here the approach of 0, assuming that it 
is always possible to improve any classical model by turning to a mixed description where 
some terms are computed using quantum mechanics. The main defect of model ^ is that 
it very often predicts too large a capture probability. To cure it, we were forced in works 
0, P] to artificially reduce the capture probability by arbitrarily reducing the capture re- 
gion. We show here that similar results can be achieved in a more self-consistent way if 
the potential-barrier-crossing process by the electron is described as a quantum mechanical 
process. In particular, two typical QM features are taken into account: (i) the fraction of 
electrons crossing the barrier from the target atom to the projectile, which classically is a 
term of exclusively geometrical origin, must be corrected by a factor ft < 1, accounting for 
the fact that a flux of quantal objects impinging on a potential hill suffers partial reflection 
even if their kinetic energy is larger than the hill's height, (ii) Furthermore, within the 
classical picture the flux of electrons to the projectile is a continuos stream while, quantum- 
mechanically, it is a resonant process, occurring only when the conditions are satisfied for 



which the electron binds to a quantized energy leveL We try to implement this feature by 
adding a modulation term, w, to the capture probability, near zero far from the resonance 
condition. 

It is important to notice that, although QM corrections are empirically added to find con- 
vergence with experiments and/or other computations, no fitting parameters are added: in- 
stead, any new parameter needed is estimated on the basis of (semi) quantitative reasoning: 
once we accept the classical-quantal mixed description, the model is entirely self-consistent. 
Another important point to stress is that our goal is to merge the QM treatment within the 
classical one without burdening too much the resulting computations. We shall show that 
through a drastical simplification of the QM computations we are able to have at the same 
time rather accurate results written in terms of simple formulas. 

Finally, we try to implement in a consistent fashion within the model the possibility by the 
target of re-capturing the electron once it has been bound to the projectile. The importance 
of this effect on the effective capture efficiency, and the limitations of the approach adopted, 
will be briefiy discussed in the appendix. 

II. DESCRIPTION OF THE OEM 

We consider a scattering experiment between a nucleus T with an active electron e, and a 
projectile nucleus P. Let r be the electron position relative to T and R the relative distance 
between T and P. Let us further consider all the three particles lying within the same plane 
V. We label the direction along the internuclear axis as the z axis and describe the position 
of the electron using cylindrical coordinates {p,z,(f) = 0). The two nuclei are considered as 
approaching at a velocity very small if compared to the orbital electron velocity. The total 
energy of the electron is 

where Zp and Zt are the effective charge of the projectile and of the target seen by the 
electron, respectively (we are considering hydrogenlike approximations for both the target 
and the projectile). Atomic units are used unless otherwise stated. Notice that our reference 
frame is non inertial, so other terms as Coriolis and centrifugal force should arise. We discard 
them on the basis of the low-velocity approximation. 



As long as the electron is bound to T we can approximate E by 



E{R) = -En-^ (2) 

with En > the unperturbed binding energy of the electron to T. 
On the plane V we can draw a section of the equipotential surface 

U{z,p,R) = E{R) = -En~^ . (3) 

The set of points (p, z) that satisfy this equation mark the limits of the region classically 
allowed to the electron. When R ^ oo this region is disconnected into two circles centered 
around each of the two nuclei. As R diminishes the two regions can eventually merge. It is 
the opening of the equipotential curve between T and P which leads to a leakage of electrons 
from one nucleus to another, and therefore to charge exchange. It is possible to solve Eq. 
(H) in the limit of vanishing width of the opening (p = 0), and find: 



Zt + yZp) — Zp 



Rm= '' ' l"" '- . (4) 
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In the region of the opening the potential has a saddle structure: along the internuclear axis 
it has a maximum at 

\/Z' 

Z = Zn = R -;= -=^ (5) 

In our work, following ||^, we assume that the electron is in a low angular momentum 
state (e.g. an s states), thus unperturbed classical electron trajectories can be visualized as 
straight segments in the radial direction, starting from the target nucleus. We assume also 
that the qualitative shape of the trajectory is not changed by the collision: the electron free 
falls towards P. 

Charge loss occurs provided that the electron is able to cross the potential barrier. Let W 
be the probability for the electron to be still bound to T at time t. We write its rate of 
change as 

^^-4/.no (a, 

Notice that we have defined the probability for an electron to be bound at P as 1 — W(t\ 
thus ruling out the possibility of ionization, which is, however, small for low-energy colli- 
sions. 



Let us now explain the meaning of the terms in the right hand side of this equation: we 
identify A^^ as the fraction of classical particles which at time t leave the target as a conse- 
quence of their motion. It is simply a geometrical factor: assuming a uniform distribution of 
electron trajectories, it is equal to the ratio between the solid angle intercepting the opening 
and the total An solid angle. The azimuthal integration is straightforward and thus 

with pm half-length of the opening in the radial direction, root of 

7-1 P t ^p / \ 

A useful approximation is to expand Eq. (||) and then Eq. (^ in powers of Pm/ R and retain 
only terms up to the second order. This was justified in the original paper ^j by the need of 
accurately modelling far collisions. However, this approximation turned out to be a rather 
accurate one since close encounters are weighted in the final cross section by the (small) 
impact parameter, thus even a rough estimate is not so much important. 
With this approximation, the calculation is straightforward and we find 



„ _ 1 ^zjz, 



( y Zp + y ZA — Zp — EnR 



(9) 



The parameter Tj is the classical period of the electron bound to T. It accounts for the fact 
that if the classical phases within the ensemble of the electrons are randomly distributed, 
during the time interval dt only the fraction dt/Tf come through the opening. The period 
can be calculated using the semiclassical relation Tt = 27™^^^ (see e.g. 0) with rics an 
effective quantum number for the electron bound to T, which can be computed from eq. (0) 
through En = Z^/{2nl^). The result is 

/ y2 \ 3/2 -, 

^ "' (i+&) 

The factor T^ = 2-n{Z'^ / {2En)Y^'^ is the unperturbed period. It is clear from eq. (p!0| ) that 
Tt/T^ < 1, and it was shown in that the average (or effective) value of this ratio is 
typically of the order of 0.5. 
Finally, ft accounts for quantum mechanical corrections to the barrier crossing probability. 
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FIG. 1: Pictorial view of the potential along the internuclear axis as it is (left figure) and how it 
is approximated by the rectangular square approximation (right figure). 

Classically ft = 1, while a flux of quantum mechanical particles impinging on a potential 
barrier is reduced by a factor /i < 1 even though the kinetic energy of the particles is larger 
than the height of the hill. In order to reduce computation of ft to an easily manageable 
form we must replace the true potential profile with a model barrier whose transmission 
factor can be analytically computed. The choice of the model potential can be crucial for 
final results, so we have tested two candidates. The former choice is the rectangular barrier; 
in fig. (|l|) we draw a schematic picture of the potential profile along the internuclear axis: 
as it is (left-hand figure), and as it is approximated (right-hand figure). The horizontal line 
labelled by "e" marks the energy of the electron. The zero of the potential well associated 
at the target is chosen so that the binding energy of the electron is equal to its unpertubed 
value En, and analogously E'^, is the binding energy to the projectile (till now undefined). 
The potential barrier between the two nuclei is depressed so that the electron is able to cross 
it with a kinetic energy AE. The width of the potential barrier is set to L, yet undefined 
but of the order of the internuclear distance R. The transmission factor TF for an electron 
coming from the left of the potential hill is 



TF 



k [{q + s) + e^'^'i {q - s)] + q [{q + s) - e^'^'^ {q - s) 



(11) 



with k = \/2En-, q = V2AE , s = ^j2E'^i electron momenta respectively in the T potential 
well, in the potential barrier, and in the P potential barrier. Of course, the relation holds 



ft 



ITFI 



(12) 



The binding energy to the projectile can be calculated by considering that, when e is bound 
to P, its energy is E{K) = ~E'^, — ZtjR, with E'^, = Zp/{n')^. At the capture radius this 
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expression and that given by eq. (^ must be equal, thus 

E'n'=En + ^^ . (13) 

We compute q along the internuclear axis: by using eqns. (|],^,|D, it is quite easy to work 
out 



It remains to estimate L. Of course, given the way eq. ( [Tl| ) has been derived, only a semi- 
quantitative estimate is required. By a straightforward application of the virial theorem, 
one finds 

and 

(y)-^*(;^) = -^n (16) 

with 2 {St) average width of the potential well. A similar relation holds for the potential well 
centered around the projectile. From the two previous equations we find {St) = l/2{Zt/En), 

{Sp) = l/2iZ,/E'J, thus, 

L = R-i<St> + <S,>)=R-^^-^^ . (17) 

Of course, we set L = when the right-hand side of the equation above is lesser than zero. 
By taking a glance at eqns. (|Tl|,|l^), one can guess that the effective number of captures 
frpNj' is strongly suppressed already for R < i?^: this is exactly what found in CTMC 
simulations (see |^). In order to take an insight at what the transmission factor looks like, 
we plot in fig. (^ ft given by Eq. (|TT|) for H - 11+ scattering. In the same way as the laws 
of quantum mechanics prevent a fraction of electron to be captured even when it would be 
classically allowed, they also-through tunnelling-would make it possible for some electrons 
to be captured even at internuclear distances R > Rm- However, it is easy to show that this 
correction to the total capture probability is very small, and thus we will neglect it. 
This model potential is extremely useful for calculations but one could wonder if it is too 
drastic an approximation, particularly in view of the fact that it is not smooth. As a test, 
we have replaced the rectangular barrier with an Eckart potential [^] 

Vix) = AVo—^^ , (18) 
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FIG. 2: Transmission factor ft (Eq. |ll|) for H - H^ collisions. Notice that the probability of 
transmission falls well below one already for values of R rather far from the maximum allowed 
capture distance Rm (Eq. Q). 

which is instead a smooth, bell-shaped curve (see fig ^. The transmission factor for this 
potential is analytically computable and, for a particle of momentum p: 

cosh(47rpA) — 1 



ft 



(19) 



cosh(47rpA) + cosh(47rv/2A2Vo - 1/16) 
The parameters \,Vo,p can be straightforwardly related to the parameters L,k,q,s given 
above. The quantity ft in (|TP|) is a mildly varying function of A but is, unfortunately, a 
strong function of E/Vq, going steeply from zero to one as this ratio crosses the unity. This 
is a deprecable feature since, as it appears clear from the equations (13-l3) above, we can 
give only semiquantitative estimates-therefore likely to have a large variability-for all the 
quantities involved. 



A main difference between the classical picture and the quantum-mechanical one is that 
the former depicts the capture process as a continuous flow. On the contrary, quantization 
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FIG. 3: Eckart potential (Eq. |l|). 

rules forbid the electronic flow from one nucleus to the other unless some resonance conditions 
are satisfied: by using the relation (|1^), we obtain 

R{n') 



Zp — Zt 
iJl E 

2 (n')2 -^" 



(20) 



that is, captures should be allowed only close to these values of i?, in correspondence of 

integer values n' . 

In this work we want to implement an algorithm to take into account the damping of capture 

probability far from these resonances, persuaded that this should drastically improve the 

predictions of the model. 

We choose to implement phenomenologically this feature by modulating ft with a weight 

function w centered around the values R{n'). For w we choose a sum of Gaussians 



w{R) 



E 



exp 



/ R - R{n') 
\AR{n')/2 



(21) 



It is necessary to find a reasonable value for AR{n'). For this we resort to the indeterminacy 
relations: infact, AR -^ means that the energetic levels for the projectile are sharply 



defined, while any finite value for AR means that they are defined only within an energy 
range AE. We suppose that the usual indeterminancy relations hold: AE x At = 1/2. 
Within the straight-line trajectory approximation, At ^ AR/u while, using again Eq. ( p!3[ ) 
and differentiating with respect to i?, we get AE k, [Zp — Zt) / R^ AR. Collecting the above 
expressions. 



The above expression breaks down when Zp = Zt. In that case we will not adopt this 
approach, instead we wil assume uniform probability of capture all along the internuclear 
distance: w = 1. Notice that, with the above definitions, it seems that we get as a bonus 
also partial probabilities for capture into well defined quantum numbers: if we are interested 
in captures into state m, it is enough to truncate the sum w (Eq. ^1]) to the single term 
corresponding to n' = m. As we shall see later, this picture however does not hold. 
Eq. (W) can be formally integrated till the end of collision: 



W{t = oo) = exp 



-^*^M^t 



(23) 



Tt 

The capture probability is P = 1 — W{oo). We assume a straight-line trajectory for the 
projectile: R = ^/W+{uty, with b impact parameter and u its velocity. Total charge 
exchange cross section is thus 

(T= 27r bPdb . (24) 

Equation (|23|) cannot be analitically integrated unless we make further simplifications. How- 
ever, it is quite easily numerically integrated by any standard mathematical software package. 

III. RESULTS 



We will benchmark the model against the experimental results from ref. |jTO[ and the 
theoretical ones coming from the molecular approach simulation of ref. [|ri[|. In fig. (^ 
we show the results for impact between multicharged hydrogenlike ions and ground state 
hydrogen. In all case impact velocity is about 1/2 (it is rigorously so for numerical results, 
while in experiments 0.49 < u < 0.51). Meyer et al. |10| published also results from other 
non-hydrogenlike ions, with the same charge states; since, for a given charge, cross sections 

10 
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FIG. 4: Charge exchange cross section versus projectile charge. Open circles, data from ref. |11]; 



open squares, data from ref. ||10|]; full squares, OBM results using rectangular barrier potential; 
triangles, OBM results using Eckart potential (Eq. O). 



are not qualitatively different, we have limited to the present selection and got rid of the 
complications involved with the definition of effective charges. 

Let us now discuss the results of fig. (^. The OBM with the choice of the rectangular 
barrier as model potential (black squares) fairly well overlaps experimental data for all cases 
but N"^"^. The reason lies in a combination of effects, which we will discuss below. Let us 
comment instead, the results from the other simulation (black triangles): they have been 
computed using Eckart potential (Eq. |T8|). Even if at first sight these results look pretty 
accurate, things are not that good. Actually, in order to reach this accuracy, we have had 
to increase the barrier height Vq so as to allow only grazing incidence, E = Vq (or g ~ 0, 
using notations of rectangular barrier potential). This, since Eckart potential has a strong 
transmission factor for E/Vq even slightly larger than one. Using the same conditions as for 
the rectangular barrier would yield an overestimate of the true cross section. 
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FIG. 5: Charge exchange cross section versus projectile charge,! < Z < 5. Open circles, data 
from ref. |11]; full squares, OBM results using rectangular barrier potential; open squares, same 
as full squares but without modulation of the transmission factor. 

Let us now consider fig. ^ it is alike fig. ^ but with 1 < Z < 5. Only data from ref. 



| TT| are available for these ions. In this case we have computes a using only ft from Eq. 
(P^) but with w given by Eq. ( PT]) (solid circles) or with w = 1 (empty circles). The only 
exception being hydrogen, for which only w = 1 can be used. It is apparent that some sort 
of modulation of the transmission factor is needed in order to avoid an overestimate of a. 
The effect of the modulation is to underestimate the cross section, but some charge states 
are more damped than others. It is the same effect, but enhanced, seen in fig. §. The reason 
of the failure of the model in these cases lies in a balance between attenuation factors and 
resonances positions. In order to better visualize it, let us consider fig. ^: there, we have 
plotted, for ions with charge Q < Z < 9, the position of each resonance R{n) [normalized 
to the maximum capture distance Rm given by Eq. (^] versus the quantum number n. 
The bars superposed to the points mark the width of the resonance: they are the Ai?(n)'s 
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FIG. 6: Normalized capture distance versus quantum number, for some highly-charged ions. 
Horizontal lines mark, for each ion the distance at which the attenuation factor ( [l^ ) is reduced to 
one half. 



given by Eq. ([2^). The horizontal lines label the distance at which the capture probability 
is reduced to 1/2 because of the factor (|T2p. Thus, capture is effective only into quantum 
numbers whose resonance position lies below this line. We compare, for example, carbon (Z 
= 6) to nitrogen (Z = 7): in both cases, attenuation ft limits capture effectively to the first 
three states, and impact parameter considerations lead to suggest that ra = 3 is the dominant 
state. However, the resonance width of n = 3 is lesser for nitrogen than for carbon, thus 
explaining the reduced capture cross section. Similar considerations, but a lesser extent, hold 
for Z = 8 and Z = 9. It is clear that slightly charged ions, for which captures only in few 
states are possible, are particularly sensitive to this effect. Different is the case for hydrogen, 
where, as explained in the previous section, no considerations of resonances apply. Notice, 
however, that oscillations in capture cross sections with the charge, due to discreteness of 
quantum levels, are not an artifact of the model but, albeit with different features, have 
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been experimentally observed |^. 

IV. CONCLUSIONS 

In this work we have presented a version of Over Barrier Model for computation of charge 
exchange cross sections. The model presents some unusual features, such as the inclusion 
within its classical background of quantum elements. We showed it to be quite reliable in a 
well defined range of parameters, namely in the low-energy, medium-to-high-projectile charge 
region, where it presents noticeable improvements with respect to other similar algorithms. 
The reasons for poorer performances in other parameters space regions have been clearly 
identified and some insights about possible further improvements are possible: the first 
candidate to work on to reach an overall fairly good accordance with experiment is the 
modulation factor w (Eq. ^1]); presently, however, we are not able to guess if a general 
optimal expression is possible, valid for all projectile-target combinations. 

APPENDIX A 

It is straightforward to adapt Eq. (j^) to accomodate the possibility for the captured 
electron to return to the target nucleus: it is enough to add a term 

^ = -Nt^JtW{t) + iV,l/,(l - W{t)) (Al) 

The second term in the r.h.s. represents the flux of electrons which have previously been 
captured by P and that now cross a second time the barrier in the opposite direction, thus 
being re-captured by T. The definition of the parameters Np, Tp, fp is the same as given in 
the previous sections, see Eqns . (|^,|TD|,|T2|) , with the trivial exchange of the projectile with 



the target. 

Eq. ( |Al|) can be solved exactly in terms of quadratures. Here, we will show that it is 

possible from this equation to recover the result for symmetrical scattering: infact, in this 

case, ft = fp, Nt = Np and Tt = Tp and, by setting ftNt/Tt = fpNp/Tp = ip for brevity, and 
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iL'fi^)^^ = "^W' we get 



1- 
1- 
1- 

l[l-e- 
2 L 



-2<I>{oo) 



-2$(oo) 



-2<I>(oo) 



1 + 



1 + 



ifir) e2*(^) dr 



! + .( 



d 

dr 



,2$(oo) 



e2*(-) dr 



1) 



(A2) 



-2$(oo) 



and we recover the prescription for equal-charge scattering. 

Of course, we can also recover easily the opposite limit, in which Zp >> Zf] in this case we 
can simply neglect the return term since fpNp/Tp « ftNt/Tf. 

We have verified that the probability of recapture is negligible for all multicharged projectiles. 
The only exception is the Zp = Zt scattering. However, we found that the use of the 
symmetrized capture probability (^) does not improve results in this case, and actually 
degrades the performances of the model. The reason is that we are overestimating the 
probability of recapture in Eq. (^): it holds rigorously in the limit of asymptotically slow 
nuclear motion, m — ^ 0; for finite m, retardation effects should be implemented to take into 
account that electrons released by the projectile at time t were captured at an earlier time 
t' . However, in the region of velocities where retardation effects can be neglected, other 
essential features of the model, such as the straight trajectory approximation, break down. 
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